Magnetic fluctuations in frustrated Laves hydrides R^Mni-^ALj^H 
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By neutron scattering, we have studied the spin correlations and spin fluctuations in frustrated 
Laves hydrides, where magnetic disorder sets in the topologically frustrated Mn lattice. Below 
the transition towards short range magnetic order, static spin clusters coexist with fluctuating and 
alsmost uncorrelated spins. The magnetic response shows a complexe lineshape, connected with the 
presence of the magnetic inhomogeneities. Rs analysis shows the existence of two different processes, 
relaxation and local excitations, for the spin fluctuations below the transition. 



I. INTRODUCTION 

There is now a growing interest in studying magnetic 
fluctuations in the paramagnetic phase of strongly corre- 
lated electron systems. An original behavior in the dy- 
namical susceptibility, where the magnetization density 
does not follow the single exponential decay predicted by 
Landau-Fermi theory is often observed. In itinerant mag- 
nets, non Fermi-liquid behavior is often associated with 
temperature anomalies of the thermodynamical proper- 
ties, like specific heat and resistivity. In Heavy-Fermions 
and superconductors, it could be the precursor a quan- 
tum phase transition around zero temperature 1 . 

In paramagnets with localized spins, a non Fermi liq- 
uid behavior of the dynamical susceptibility could be also 
expected due to strong magnetic disorder, since magnetic 
inhomogeneities associated with disorder and frustration 
yield a distribution of relaxation times. However, in usual 
spin glasses with localized spins, although such distribu- 
tion indeed exists, the excitation spectrum in the para- 
magnetic phase usually remains much narrower than in 
itinerant magnets, yielding typical energies of the fluctu- 
ating spins below 1 meV. 

The intermctallic Laves compounds RMn2 have two 
main characteristics 2-4 : i) they are at the frontier be- 
tween localized and itinerant magnetism, the instability 
of the Mn spins being governed by the distance between 
first neighbor Mn pairs, and ii) they show a topologi- 
cal frustration of the pyrochlore Mn lattice for antifer- 
romagnetic first neighbor Mn-Mn interactions. When 
R is a magnetic rare earth, the diamond-like R lattice 
interacts with the frustrated Mn lattice. Varying Mn- 
Mn interatomic distance d allows one to encompass the 
critical distance do between localized and itinerant Mn 
magnetism. For d > do, spontaneous Mn moments are 
present and the frustrated Mn-Mn interactions dominate, 
whereas for d < do Mn moments are induced by rare 
earth moments and the topological frustration is sup- 
pressed. Changing interatomic distances is realized by 
chemical substitution 2-4 or by applying pressure 5 . 

Introducing hydrogen in interstitial sites also modifies 
the magnetic properties 6 . In Laves compounds, hydrogen 



expands the lattice, favoring the localization of the Mn 
moments, but also modifies the frustration by making 
Mn-Mn interactions nonequivalent. Therefore, ordered 
H atoms release the topological frustration, yielding long 
range ordered magnetic structures 7-10 . On the other 
hand, H disorder induces a random variation of the ex- 
change interactions and therefore magnetic disorder 1 ' 12 . 

In Y(Mni„ a; Al a; )2Hy and Y(Mni_ x Fe a; )2H y , we stud- 
ied such phases by magnetic neutron diffraction 11,12 . 
Surprisingly, and in contrast with usual spin glasses, we 
found that these disordered magnetic phases combine 
short spin correlation lengths (around 10-20 A) with high 
transition temperatures (100-250 K), in the same range 
as in the ordered parent compound. The behavior of 
the short range order parameter was intermediate be- 
tween that expected for a canonical Edwards-Anderson 
spin glass, and the Brillouin curve for a non frustrated 
magnet. 

Up to now the study of the spin fluctuations and 
magnetic transition in such topologically disordered 
phases has received little attention. A previous study 
of Y(Mno.gAlo.i)2 single crystal by inelastic neutron 
scattering 13 , where interatomic Mn distances are close 
to the instability threshold, revealed very interesting fea- 
tures with regards to canonical spin glasses. Namely, it 
pointed out the existence of a large energy range of the 
magnetic fluctuations (up to about 40 meV), well above 
the typical fluctuation energy in spin glasses, and close 
to that observed in other itinerant magnets 14 ' 15 . This 
large energy range was attributed to the influence of both 
topological frustration and itinerant character of the Mn 
moments, without any way to clear up their respective 
roles. 

Studying the magnetic excitations yields precious in- 
formation about the microscopic nature of the magnetic 
ground state and ordering transition. Especially, it could 
help to explain why the transition temperature remains 
high in spite of the strong disorder. Inelastic neutron 
scattering is the most straightforward way to probe such 
excitations. We have studied the spin fluctuations in 
R(Mni_ 2; Al 2; )rIy compounds (R = Y or Tb) by time of 
flight neutron scattering, together the steady-state cor- 
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relations by neutron diffraction, and the low field sus- 
ceptibility by static magnetization measurements. In all 
compounds, H disorder was achieved by substituting a 
small amount of Al on the Mn sites. By choosing appro- 
priate compounds, we checked the influence of magnetic 
moments on the R sites, hydrogen concentration and hy- 
drogen disorder. The results suggest an original picture 
of a 'cluster glass' magnetic state. 

II. SAMPLE DESCRIPTION 

Powdered samples of R(Mni_ x Al x )2Dj / , with (R = Y, 
Tb), Al content £=0.09, and different deuterium con- 
tents were prepared using the technique described in Rcf. 
11, and characterized by X-ray diffraction. All samples 
showed a single chemical phase, besides a very small 
amount of RH3 impurity phase in some cases. We mostly 
studied three disordered samples, Tb(Mn .gAlo.i)2Dj,, 
(y=l and j/=4), and Y(Mn . 9 Al .i)2D !y (y=4). Exact D 
contents are given in Table I. In most cases, the deu- 
terium isotope was chosen instead of hydrogen to de- 
crease the incoherent background. The small Al con- 
tent induces hydrogen disorder and breaks down the long 
range magnetic order 11 . The Tb ion with high magnetic 
moment yields a strong enhancement of the magnetic 
fluctuations (by about a factor 10), yielding an excellent 
precision on the lineshape of the dynamical susceptibil- 
ity. The influence of this second magnetic ion on the 
spin correlations and spin dynamics was checked in com- 
parison with Y(Mno.gAlo. 1)2^4 where Y is non magnetic. 
The amount of hydrogen (or deuterium) introduced (y=l 
or y=4) allows to vary interatomic distances in a large 
range. In this range, the djyrn-Mn distance (2.77 A for 
Tb(Mn . 9 Alo.i) 2 Di, and 2.87 A for TtyMno.gAlo.i^D^ 
remains above the magnetic instability threshold for the 
Mn moments (do— 2.66 A). 

In comparison with the disordered compounds we 
studied the spin fluctuations in ordered TbMn 2 D 4 , 
TbMn2H4, and YM112H4. Their simple long range an- 
tiferromagnctic structure, which is determined by the H 
superstructure 7,8 , is a reference to check the influence of 
hydrogen or deuterium disorder on the spin fluctuations. 



III. NEUTRON DIFFRACTIION 
MEASUREMENTS 

We studied the magnetic correlations on the powder 
diffractometer G6-1 of the reactor Orphee (with an inci- 
dent neutron wavelength of 4.734 A). Temperature var- 
ied between 350 K and 10 K. A spectrum measured in 
the paramagnetic range was subtracted to extract the 
magnetic correlations. The magnetic spectra of the Al 
substituted samples show diffuse magnetic peaks, corre- 
sponding to short range antiferromagnetic correlations 
(Fig. 1). The peak positions are indexed with the same 



propagation vector as in the ordered parent compound 
RMn2D4. Using the Fullprof program 16 , wc refined the 
data assuming short range ordered antiferromagnetic re- 
gions of diameter L c with ordered moments /zr and /iMn, 
keeping the same type of order as in Ref. 8. Results are 
given in Table I. We find correlated regions of exten- 
sion L c between 15 and 30 A. L c is temperature inde- 
pendent, but varies with H content and R substitution. 
The presence of magnetic Tb ions increases L c with re- 
spect to non magnetic Y. It also induces ferromagnetic 
correlations on a shorter length scale (about 6 A). These 
ferromagnetic correlations, shown by an increase of the 
magnetic scattering at low angles coexist with the anti- 
ferromagnetic ones. The short range ordered moments 
remain well below the moments in the long range or- 
dered samples (/^Mn=4 [Ib and /ZTt>=9 [Ib in ordered 
TbMn2D4 for instance). In spite of the strong magnetic 
disorder, the transition temperatures Tsro towards the 
short range ordered state are in the range 100-250 K, 
therefore not much decreased with respect to the Neel 
temperature of the long range ordered state (Tn =300 K 
in TbMn 2 D 4 and 330 K in YMn 2 D 4 ). This was also the 
case of other disordered hydrides previously studied 12 . A 
detailed description of the short range magnetic order in 
the R(Mn .9Al .i)2H y family will be given later. 

IV. MAGNETIC MEASUREMENTS 

The magnetization was measured with a SQUID mag- 
netometer, under static magnetic field of 10 Oe, in the 
zero field cooling (zfc) and field cooling (fc) conditions. 
Typical curves are shown in Fig. 2. In all disordered 
samples, magnetic irreversibilities start around 230-280 
K, namely above Tsro- In the Y-samples, where only 
Mn sublatticc is magnetic, the fc curve continuously in- 
creases with decreasing temperature, well above the zfc 
curve, which shows a broad maximum then flattens or 
decreases. The magnetization of the two Y-samples have 
the same order of magnitude in the paramagnetic regime, 
but in the irreversibility region, the magnetization is 
much smaller for the sample with low hydrogen content. 
The temperature dependence of both fc and zfc curves is 
reminiscent of the behavior in cluster glasses 17-19 , where 
magnetic order occurs within spatially limited regions, 
and the relative orientations of the superparamagnctic- 
like spin clusters depend on the cluster interaction field, 
anisotropy field, and thermal activation 20 . In the Tb- 
samples with two magnetic sublattices, this cluster glass 
behavior is superimposed on a global paramagnetic com- 
ponent, with a variation close to 1/T. This second com- 
ponent should be related with the strong thermal fluc- 
tuations of the Tb moments, which become progres- 
sively polarized by the Tb-Mn exchange field as temper- 
ature decreases 8 . In all the short range ordered samples, 
anomalies are observed in the low temperature range (30- 
50 K), as shown by a small maximum in the zfc curve, 
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and an upturn of the fc curve. Similar anomalies were 
predicted and observed in some cluster glasses, when the 
local anisotropy fields start to be significant 20 . They 
suggest the presence of small superparamagnetic clus- 
ters which could fluctuate down to low temperatures, 
having their anisotropy fields perpendicular to the in- 
teraction fields 20 . In Y(Mn .gAlo.i)2D4, the low temper- 
ature anomalies result in a plateau on the fc curve. This 
plateau, whose origin is unexplained, reaches only about 
3 % of the maximum expected due to demagnetization 
effect. 

In comparison with the short range ordered com- 
pounds, the magnetization of YMJI2D4 and TbMn2D4 
with long range antiferromagnetic order is shown in the 
insets. In both samples, the magnetic irreversibilities 
start abruptly at the Neel temperature Tn measured by 
neutrons 7,8 . In YM112D4, a magnetization plateau is ob- 
served below Tn, with a value of about 3 % of the de- 
magnetization value. This plateau, rather unusual for 
an antiferromagnet, is presumably connected with the 
first order character of the transition, yielding a spin 
blocking induced by hydrogen order and ordered anti- 
ferromagnetic moment which remains close to saturation 
almost up to Tn- TbMn2D4 shows a similar behavior, 
with an additional paramagnetic component. In both 
samples, the low temperature anomalies observed in the 
short range ordered samples are absent. 



V. INELASTIC NEUTRON SCATTERING 

A. experimental set up and data analysis 

The experiments were performed in the time of flight 
spectrometer Mibemol of the reactor Orphee, with an in- 
cident wavelength of 5 A (incident energy of 3.27 meV). 
The energy resolution (FWHM=0.15 meV) and the effi- 
ciency of the detectors were determined by measuring a 
vanadium standard of the same plate-like shape as the 
samples. After background and absorption corrections, 
the neutron cross sections were calibrated in absolute 
scale using the vanadium spectrum. Measurements were 
performed versus temperature, using cither a closed cycle 
nitrogen flow cryo-oven (between 350 K and 100 K) or 
a helium flow cryostat (between 300 K and 20 K). The 
maximum temperature of 350 K was chosen to avoid any 
decomposition of the samples due to hydrogen migration. 
The temperature range allowed us to study the spin fluc- 
tuations in the paramagnetic range, and to follow their 
evolution below the transition. The range of the elastic 
scattering vector is the same as in the diffraction mea- 
surements. 

The scattering cross section for magnetic quasielastic 
neutron scattering can be expressed in a general form 



S(q,u ) ) = -n(u) X '^) 

7T 



(1) 



where n(u>) is the Bose factor taking into account the 
thermal population factor of the fluctuations, and Xqi 1 ^) 
the dynamical susceptibility yielding the response of the 
system. Within linear response theory, \" q ( w ) 1S the prod- 
uct of the static transverse susceptibility x(q,w=0), and 
the spectral function F(uj). In the simplest case of a sin- 
gle exponential decay (SED), there is a single energy T 
(or equivalently a single relaxation rate T=l/27iT, taking 
h=l) for the magnetic fluctuations, so that the spectral 
function is expressed as F(w) = ^ fsi^p. ■ In classical spin 
glasses like CuMn, the spectral response was expressed 
by a distribution of Lorcntzians correspondig to a dis- 
tribution of relaxation rates 23 ' 24 . Still, the quasielatic 
scattering cross section in the experimental energy range 
was well fitted by a single Lorentzian lineshape 23 . 

This simple model could not account for our data, 
yielding strong misfits even in the paramagnetic range. 
Therefore we considered explicitely an energy distribu- 
tion in the spectral lineshape. We took two limiting 
forms for this distribution, namely a bimodal distribution 
(model 1) and a uniform distribution (model 2). Model 1 
assumes two typical energies Ti and T 2 (ri<r 2 ) for the 
spin fluctuations. It was used previously to analyze the 
spin dynamics in cluster glasses 21 ' 22 . Model 2 assumes 
that the energies of the spin fluctuations are distributed 
with an equal probability between lower and higher en- 
ergy bounds Ti and T 2 . It was recently proposed by 
Bernhoeft 25 within a general model of the dynamical 
phase inhomogeneities in non Fermi liquid systems. 

The quasielastic neutron cross section is written 

- within model 1: 



S{q,u) = 



j/kT _ 1 



(Fi(w)Xi(g) + F 2 (w)x 2 («)) (2) 



where xi an d X2 are the static susceptibilities of the two 
quasielastic contributions and the spectral functions are 
respectively 



1 ri 



1 



(3) 



with Ti, T2 being the corresponding energy lincwidths. 
The static susceptibility is defined as x= Xi+X2- 
- within model 2: 



(4) 



where x is the static susceptibility and F(u>) the spectral 
function written as 



I ht^-ht,) I ( arctan (rT) - -ctan(^) 



(5) 



where Ti and T 2 are respectively the low and high energy 
limits of the energy distribution. 

In the above expressions, w>0 corresponds to a neu- 
tron energy gain. The prefactors in equations (3) and 
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(5) normalize to unity the energy integration of the spec- 
tral function. Therefore, in the high temperature limit 
w-CkT (valid here in the paramagnetic range and down 
to about 50 K), the quantity T% is equal to the energy- 
integrated intensity of the quasielastic signal. The fitted 
energy widths do not depend on the normalisation factor. 

In the limit of a Curie- Weiss behavior, the neutron sus- 
ceptibility is related to the effective fluctuating moment 

2 2 

by the expression: x(q=0,T)=A^^-^-. A corresponds 
to the calibration of the neutron cross section in abso- 
lute scale. This expression provides an effective moment 
value, to be compared with the value from bulk magnetic 
measurement. 

The total cross section is written as 

S tota \q,uj)=c el S(Lu) + S(q,u J ) (6) 

The neutron intensity I detected versus time of flight t 
at the scattering angle 29 is proportional to the quantity: 

/(M) « is tota W) (7) 

The neutron intensity was convoluted with the exper- 
imental resolution and fitted to each individual time of 
flight (TOF) spectrum, with c e i (intensity of the elastic 
peak), r 1; T 2 , and xi, X2 (model 1) or x (model 2) as pa- 
rameters. A very small linear background was also added 
in the fit. The fit was performed in the energy range (- 
20, 2) meV, where the phonon contribution is small and 
could be neglected. 

B. results 

a) short range ordered samples 

We first show the results in Tb(Mn .9Al .i) 2 D y (y=l 
and y=4). In spite of their different D contents, ordering 
temperatures and magnetic correlations lengths, the be- 
havior of the magnetic fluctuations is the same for the two 
samples. Typical TOF spectra are shown in Figs. 3 and 
4 for temperatures respectively above and below Tsro- 
Note that due to the t~ 4 factor of equation (7), the max- 
ima of the quasielastic signals are shifted with respect 
to the elastic one. In spite of this distortion, the time 
dependent plot were preferred to the energy dependent 
ones (in inset), since they provide a better estimation of 
the fit quality at low energies. 

Whatever the model chosen for the spectral function, 
the general evolution with temperature can be described 
as follows. Above Tsro, the susceptibility follows a Curie 
law, and the quasielastic integrated intensity T% is tem- 
perature independent (Fig. 5). Below Tsro, we observe 
a gradual transfer of intensity from the quasielastic sig- 
nal to the elastic one, the sum of the two contributions 
being conserved with temperature. The high energy sig- 
nal of width T 2 dominates, with a weight of 80% of the 
succptibility at high temperature. The energy linewidths 



have a rather original behavior (Fig. 6). Above Tsro, IT 
and r 2 are temperature independent. Below Tsro, the 
small linewidth r 4 starts to decrease whereas the large 
linewidth T 2 increases with decreasing temperature. This 
behavior is the same for all q values. In the paramagnetic 
regime, the quasielastic cross section is well fitted by the 
two models. Below Tsro, model 1 still yields very good 
fit of the data, whereas model 2 starts to show some mis- 
fits. They occur mainly below 80 K where the linewidths 
i\ and T 2 are very different. We have therefore restricted 
the use of model 2 to temperatures above 80 K, where 
this model was meaningful. In the paramagnetic range, 
the two models give an energy linewidth in the same en- 
ergy range, with Ti around 0.1 meV, and T 2 around 3-5 
mcV. 

The q dependence of the susceptibility and quasielas- 
tic widths is shown in Fig. 7. The susceptibility x shows 
a weak and broad modulation with respect to the form 
factor of the free ion. This modulation remains the same 
below Tsro- The magnetic correlations observed below 
Tsro by neutron diffraction correspond to a modulation 
of the elastic component. The energy linewidths are al- 
most q- independent, both above and below Tsro- The 
Tb(Mno.gAlo.i) 2 Di sample, where the magnetic correla- 
tions are enhanced, shows the same behavior. 

In Y(Mn .9Alo.i) 2 D 4 , the magnetic signal coming from 
Mn ions only is much smaller, but still shows the complex 
lincshape described above. The main difference with the 
Tb samples concerns the temperature dependence of the 
large energy linewidth T 2 , which now strongly decreases 
below Tsro (Fig. 8). This is observed whatever the 
model used to describe the spectral function. 

b) long range ordered samples 

In comparison, we show the magnetic fluctuations in 
TbMn 2 D 4 , TbMn 2 H 4 , and YMn 2 H 4 where long range 
magnetic order occurs at Tn through a first order magne- 
tostructural transition. Above Tn , the magnetic fluctua- 
tions are very similar in the ordered and in the disordered 
samples. The main effect of the LRO is a strong suppres- 
sion of the quasielastic intensity below Tn- For example, 
at T= Tn/2 the intensity of the magnetic fluctuations 
is decreased by a factor 3 with respect to its value in 
the paramagnetic state in the ordered Tb samples, and 
by a factor 5 in ordered YMn 2 H 4 , whereas it decreases 
by 30% in disordered Tb(Mn .gAl .i) 2 Di and only 8% in 
Y(Mno.gAlo.i) 2 D 4 . The second important difference is 
the occurence of a well defined localized excitation in the 
ordered Tb samples below Tw- 
in Fig. 9, we show two typical TOF spectra for 
TbMn 2 D 4 at 340 K and 160 K, namely above and below 
Tn- At 340 K, the spectral lineshape is well fitted by a bi- 
modal distribution (model 1), with typical energy widths 
of 0.5 and 2 meV, therefore comparable to the energy 
widths of the disordered samples at the same tempera- 
ture. Below Tn, a purely inelastic signal appears in the 
TOF spectra, coexisting with the quasielastic one. This 
new component was fitted by a Gaussian lineshape, of 
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energy Eo and width Ti nc \. Its intensity (fig. 10a) shows 
a broad maximum versus temperature. The energy Eo 
slightly increases from 6 to 7 meV when temperature de- 
creases (Fig 10b), with constant energy width rj no i of 2 
meV. The quasielastic signal strongly decreases just be- 
low Tn, but still persist in the ordered state. The large 
energy width increases with decreasing T, whereas the 
small energy width Ti is about temperature independent. 
The sum of the three contributions (elastic, quasielastic 
and inelastic) is conserved with temperature. 

In TbMn 2 H 4 , the inelastic and quasielastic signals are 
the same as in the TbMn2D4. The elastic component is 
enhanced by about a factor 10 in TbMn 2 H4, due to the 
large incoherent contribution of the hydrogen isotope. As 
expected from previous inelastic neutron studies of Laves 
hydrides 26 , the dynamics of hydrogen ordering cannot 
be resolved in the present experiment. In other words, 
hydrogen (or deuterium) is frozen at the time scale of 
the neutron experiment (t=9x 10~ 12 s), and there is no 
interference between H or D diffusion and the magnetic 
fluctuations. 

Finally, in ordered YMn 2 H 4 (Fig 11), the quasielastic 
signal above Tn also consists of two components. Below 
Tn, it becomes very small but can still be detected down 
to 100 K. The energy width in the paramagnetic state are 
about the same as in the disordered Y(Mno.gAlo. 1)204, 
but here they do not vary much with temperature. The 
inelastic localized mode observed in the ordered Tb sam- 
ples was not observed. 

VI. DISCUSSION 

A. paramagnetic regime 

In the high temperature range where the samples show 
a Curie- Weiss behavior, the effective moments deduced 
from the integration of the quasielastic neutron cross sec- 
tion agree with the value deduced from magnetization 
(table II). This means that most of the spin fluctua- 
tions are probed by the experimental time window, and 
rules out the presence of high energy spin fluctuations 
(in the 50 meV range or above), similar to that observed 
in YMn 2 27 ' 28 , Y(Mn .gAlo.i)2 13 or other itinerant sys- 
tems. Therefore, all samples are in the localized moment 
limit, even at low H(D) content, and the wide energy 
distribution seen in both Tb and Y samples should be 
connected with peculiar features of the magnetic disor- 
der. At high temperature, for localized spins coupled by 
near neighbor interactions, the energy linewidth of the 
magnetic fluctuations is related to the values of the ex- 
change constants 29 . This is the case of the present sam- 
ples, where the anisotropy energy is much smaller than 
the exchange energy 30 . Within molecular approximation, 
the lineshape of the response function is expected to be 
intermediate between Gaussian and Lorentzian, and its 



second moment, which measures the quasielastic width 
is given for a polycrystallinc sample by 

Z?= (8/3),S(^+l)^J 2 (R J )(l-sin(qR l )/(qR 4 )) (8) 

i 

In this model, the energy linewidth is temper- 

ature independent, in agreement with experiment, and 
in contrast with the behavior of itinerant systems where 
the linewidth increases linearly with temperature (Ko- 
rringa behavior). Assuming first neighbor interactions 
only (Ri=2.86 A), the smooth variation of the linewidth 
predicted by equation (8) above 1 A -1 is also consistent 
with experiment (Fig. 7b), although the expected de- 
crease at low q values is not seen. 

The distribution of linewidths may be connected with a 
distribution of exchange constants, induced by the frozen 
hydrogen or deuterium disorder. As discussed previ- 
ously for the ordered compounds YMn 2 D 4 , the pres- 
ence of a H(D) atom near a Mn-Mn pair modifies the 
amplitude, and even the sign, of the first neighbor ex- 
change interaction 8 . The various environments of the 
Mn atoms having zero or one H(D) atom nearby the 
Mn pairs should be at the origin of the linewidth dis- 
tribution. In the Tb samples, the magnetic fluctuations 
are dominated by the Tb moments of much larger mag- 
nitude. The distribution of linewidths therefore mostly 
measures the distribution of Tb-Mn interactions, showing 
that they are also sensitive to the number of H(D) near 
neighbors atoms. This is also consistent with the analysis 
of the ordered structure 8 . Taking 5mii=2, and summing 
equation (8) over the 6 first neighbors of a Mn ion, we 
obtain from the two energy widths in Y(Mn .9Al . 1)204, 
low and high limits for the Mn-Mn exchange constant 

\ J 2 Mn-Mn equal to 0.03 and 0.3 meV. A similar evalu- 
ation for Tb(Mn .gAlo. 1)204 (<Srb = 3) yields limit values 

for the the Tb-Mn exchange constant \J J 2 Tb-Mn equal 
to 0.03 and 0.1 meV. 



B. ordering transition 

In the long range ordered compounds, the transition 
is of first order, whereas in the short range ordered com- 
pounds, it is close to second order 7,11 . As shown pre- 
viously, H(D) order is the key parameter which controls 
the frustration. In the LRO samples, the first order tran- 
sition is directly connected with the formation of H(D) 
superstructure which relieves the frustration of the Mn 
lattice. In contrast, in the disordered samples, the second 
order transition occurs in a range of temperature where 
the slowing down the H(D) diffusion yields a more pro- 
gressive change from thermal to frozen disorder. In both 
cases, the present results show that the transition is char- 
acterized by a deep change in the magnetic fluctuations, 
namely, i) a transfer of intensity from the quasielastic to 
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elastic signal connected with the onset of static correla- 
tions, and ii) a change in the behavior of the quasielastic 
widths, which start to vary with temperature. 

The transition temperature Tsro remains in the range 
100-250 K, therefore is not much decreased with respect 
to the Neel temperature of the ordered compounds (300- 
330 K). Such high transition temperature is often ob- 
served in cluster glasses 31 . It contrasts with classical 
spin glasses where the spin glass temperature is typicality 
ten times smaller than the transition temperature of the 
ordered parent compound 31-33 . The average exchange 
constant Jiin-Mn deduced from Tn in the mean field ap- 
proximation (T N =X) j |JMn-Mn'S'Mn('S'Mn + 1), summing 
over the 6 first neighbors of a Mn ion) yields JMn-Mn=0.9 
meV, which is about 3 times higher than the upper limit 

for \J J 2 Mn-Mn deduced from the linewidths in the para- 
magnetic regime. Obviously, such evaluations only pro- 
vide an order of magnitude of the exchange constant since 
they neglect the influence of inhomogeneity and frusta- 
tion. 

A simple picture in terms of a cluster glass can ac- 
count for the observed decrease of Tn with respect to 
the ordered state. One assumes clusters of diameter L c 
(where L c is the correlation length) consisting of an or- 
dered core surrounded by a disordered layer of thickness 
d (where d=2.8A is the first neighbor distance). In first 
approximation, the decrease of Tn is given by the per- 
centage disordered spins, assuming that they do not con- 
tribute to the average exchange interaction. For the Tb 
samples (Tn=300K in ordered TbMn2D4), this simple 
picture yields short range ordering temperatures of 112K 
(y=4) and 161K (y=l), in relatively good agreement with 
the experimental values (135K and 200K respectively). 
A better evaluation would require to define different in- 
ter and intracluster exchange interactions, as proposed 
in mean fields models of cluster glasses 34 . 

The coexistence of highly frustrated and low frustrated 
spins yields a large distribution of energy barriers, as ob- 
served in cluster glasses or inhomogeneous superparam- 
agnets. The magnetization data support this qualita- 
tive picture, showing that some spins could freeze either 
above Tsro (up to the Neel temperature of the ordered 
compound), or well below Tsro down to the low T range. 

C. spin dynamics below the transition 

As temperature decreases, more and more spins freeze 
and contribute to the static clusters. The fluctuating 
spins remain almost uncorrelated like in the paramag- 
netic regime. In Y(Mno.gAlo. 1)204, the decrease of the 
two quasielastic widths suggests that these uncorrelated 
spins relax between several metastable configurations, 
with typical times (t=1/27iT) increasing with decreas- 
ing temperature. 

In the Tb samples, the formation of the static clusters 
is associated with the emergence of two types of spin dy- 



namics : a relaxation process involving weakly coupled 
and almost uncorrelated spins, as shown by the decrease 
of the small quasielastic width T\, and an excitation pro- 
cess as shown by the increase of the large quasielastic 
width T 2 . Here the magnetic signal is dominated by the 
Tb spins. In ordered TbMn2D4, a local mode of energy 
E is clearly observed below Tn, and is attributed to local 
fluctuations of the Tb spins in the static Mn molecular 
field. Such mode is predicted and observed in ordered 
RFe2 or RA1 2 Laves phases 35 ' 36 . Its energy is expected 
to increase smoothly with temperature like the ordered 
Mn momcnt(Eo oc Jxb-Mn Sm-o)- In the short range or- 
dered Tb compounds, one should also expect a similar 
inelastic mode, but broadened by the distribution of ex- 
change interactions and situated at a lower energy, due 
to the lower value of the short range ordered Mn mo- 
ment. Therefore, this mode could be hidden by the large 
quasielastic signal. In order to check this assumption, we 
have performed for Tb(Mno. 9 Alo. 1)21)4 the same fit as for 
ordered TbMn2D4. It yields a slight improvement of the 
fit quality, with an inelastic contribution of about 10% 
of the quasielastic one. The energy E is smaller than 
in TbMn2D4 and increases more rapidly with decreasing 
temperature (from 3 to 5 meV), as expected from a sim- 
ple model. The quasielastic widths, especially the low 
energy width Ti are not much affected. 

This analysis points out the strong similarities between 
the spin fluctuations in the short range and long range 
ordered Tb compounds. It also suggests that the large 
quasielastic signal and the inelastic one have the same 
nature, and we note that their characteristic energies Eo 
and T2 vary in the same way with temperature. The 
onset of short range order weakens the energy of the in- 
elastic mode, and changes the balance between the two 
contributions. The persistence of a quasielastic signal in 
the ordered state, which is not observed in RFC2 Laves 
phases, is presumably related to the influence of hydrogen 
or deuterium. The lowering of symmetry due to the H(D) 
superstructure 8 might create lower energy levels which 
cannot be distinguished from the quasielastic scattering. 
The residual disorder, since one H(D) site is halfly occu- 
pied, could also play a role. Whatever the exact mech- 
anism, we have observed local excitations of low energy 
which develop as soon as a static magnetic order freezes 
in, and are almost insensitive to the correlation length. A 
better analysis of these excitations would require higher 
resolution measurements on a single crystal. 

VII. CONCLUSION 

In the Laves compounds R(Mno.gAlo.i)2Dj, with deu- 
terium or hydrogen disorder, a transition towards short 
range magnetic order is probed by neutron scattering. It 
is shown by the occurrence of static clusters of a few tens 
of Angstroms, with a size independent of temperature. 
The transition temperature Tsro is high in spite of the 
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small cluster size, and the transition is close to second or- 
der, in contrast with the first order transition in the long 
range ordered compounds. In the paramagnetic range, 
the complex lineshape of the spectral function suggests 
a distribution of exchange constants, possibly related to 
the onset of noncquivalcnt magnetic bonds and the in- 
fluence of frozen H(D) configurations. Below Tsro, sev- 
eral types of fluctuations are observed, depending on the 
presence of the Tb magnetic ion. They are attributed 
either to the relaxation of almost uncorrelated spins, or 
to local spin excitations within the frozen clusters. The 
spin fluctuations in the short range and long range or- 
dered compounds show strong similarities. The overall 
behavior and the macroscopic magnetic properties are 
interpreted in terms of a cluster glass, where short range 
ordered regions with local H(D) order and low frustra- 
tion coexist with strongly frustrated and uncorrelated 
spins. It would be interesting to know if the observed 
spin fluctuations arise from the peculiar coupling of hy- 
drogen and magnetic lattices, or if they could be observed 
in other cluster glasses or disordered magnets with topo- 
logical frustration. 
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Figure captions 

Fig.l: Magnetic neutron diffraction spectra at 
10 K for Tb(Mn .9Alo.i) 2 Dj / {y=A and y=l) and 
Y(Mn .gAlo. i)2^a, measured on G6-1 diffractometer (A 
=4.734 A). A spectrum measured in the paramagnetic 
range was subtracted to extract the magnetic contribu- 
tion. For Tb(Mn .9Al .i)2D4, the first antifcrromagnetic 
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peak is not observed since magnetic moments are along 
direction <111>. Solid lines are fits of the data, yielding 
the short range ordered magnetic moments and correla- 
tions lengths given in Table 1. 

Fig. 2: Low field magnetization in Tb(Mno.gAlo.i)D 3/ 
(y=4 and y=l), and in Y(Mno.gAl .i)D !y (y=4 and y=l), 
measured in a static field H=10 Oe, using zero field cool- 
ing (zfc) and field cooling (fc) conditions. Insets: low 
field magnetization of the long range ordered YMn 2 D 4 
and TbMn2D4 compounds. The transition temperature 
Tsro towards short range order is shown by arrows. 
Other arrows point out the low temperature anomalies. 

Fig.3: 

Time of flight (TOF) spectrum of Tb(Mno. 9 Alo.i)D 4 at 
250 K (paramagnetic range), focusing on the quasielas- 
tic region. Spectra were regrouped for detectors within 
an angular range 23.5-63.1 degrees. Fits with model 1 
(a) and model 2 (b) are shown by solid lines. Inset: the 
corresponding neutron cross section versus energy trans- 
fer, and a schematic drawing of the distribution of energy 
widths for each model. 

Fig.4: TOF spectrum of Tb(Mno. 9 Alo.i)D 4 at 80 K 
(below Tsro), in an angular range 23.5-63.1 degrees. 
Fits with model 1 (a) and model 2 (b) are shown by 
solid lines. Inset: the corresponding neutron cross sec- 
tion versus energy transfer and a schematic drawing of 
the distribution of energy widths for each model. 

Fig. 5: Elastic, quasielastic intensity T%, and 
sum of the two contributions versus temperature in 
Tb(Mn .9Al .i)2D4- TOF spectra have been regrouped 
within an angular range 23.5-63.1 degrees. Models 1 and 
2 give equivalent results. 

Fig. 6: Quasielastic widths Ti and T 2 versus tempera- 
ture in Tb(Mn .gAlo.i)2D4, obtained by fitting the data 
with model 1 (a) and model 2 (b). 

Fig. 7: Results obtained in Tb(Mno.gAlo. 1)204 with 
model 1. a) Quasielastic intensity versus q for two tem- 
peratures, 200 K and 80 K, above and below Tsro re- 
spectively. Solid lines correspond to the squared mag- 
netic form factors of Tb and Mn ions. Dotted lines are 
guides to the eye. b) Quasielastic widths Ti and T 2 ver- 
sus q at 200 K. Solid lines are fits with the equation (8) 
described in text, c) Quasielastic widths versus q at 80 
K. 

Fig. 8: Quasielastic widths Ti and T 2 versus tempera- 
ture in Y(Mn .gAl . i)2D4, obtained with model 1. 

Fig.9: TbMn 2 D 4 : TOF spectra at 340 K and 160 K, 
(Tn =300 K). In inset, the corresponding neutron cross 
section versus energy transfer. Solid lines are fits with 
model 1, adding an inelastic contribution below Tn- 



Fig. 10: TbMn2D 4 : a) variation of the integrated in- 
tensities versus temperature, b) quasielastic widths Ti 
and F 2 and energy E of the inelastic signal versus tem- 
perature. 

Fig.ll: YMn2H 4 : a) variation of the integrated inten- 
sities versus temperature, b) quasielastic widths Ti and 
T 2 versus temperature. 



Table captions 



R 


y (exact value) 


Tsro (K) 


R (A) 


fJ-Mn (/Us) 


A*R (Ms) 


Tb 


4 (3.87) 


135 (10) 


20 (2) 


1.7 (3) 


3.5 (5) 


Tb 


1 (1.04) 


200 (10) 


30 (3) 


2.5 (3) 


8.8 (8) 


Y 


4 (3.81) 


250 (10) 


13 (2) 


1.8 (2) 




Y 


1 (1.38) 


190 (10) 


15 (3) 


3.0 (3) 





TABLE I. Lransition temperature Tsro, correlation 
length L c and short range ordered magnetic moments at low 
temperature (8-11 K) deduced from neutron diffraction data 
in the R(Mn .9Al .i)2Dj / samples. Error bars are mentioned 
in units of the last digit. 



Compound 


Effective magnetic moment (/is) 
(Magnetization) (Neutron scattering) 


Y(Mno. 9 Alo.i)2D 4 
Tb(Mno. 9 Alo.i)2D 4 
Tb(Mn . 9 Alo.i)2Di 


4.1 ± 0.2 4.0 ± 0.6 

10.1 ± 0.5 11.9 ± 1.8 

10.2 ± 0.5 11.6 ± 1.7 



TABLE II. Effective magnetic moments per chemical fo- 
mula in the paramagnetic regime, deduced from magnetic and 
inelastic neutron scattering measurements. 
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